Ising antiferromagnet with ultracold bosonic mixtures confined in a harmonic trap 
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We present accurate results based on Quantum Monte Carlo simulations of two-component bosonic 
systems on a square lattice and in the presence of an external harmonic confinement. Starting from 
hopping parameters and interaction strengths which stabilize the Ising antiferromagnetic phase in 
the homogeneous case and at half integer filling factor, we study how the presence of the harmonic 
confinement challenge the realization of such phase. We consider realistic trapping frequencies 
and number of particles, and establish under which conditions, i.e. total number of particles and 
population imbalance, the antiferromagnetic phase can be observed in the trap. 
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In recent years great attention has been devoted to 
mixtures of ultracold atomic gases in optical lattices [T]- 
8 . One of the most remarkable features of such systems is 
the possibility of stabilizing quantum magnetic phases [9, 
IIP] , Though interesting and fundamental on their own, 
better understanding of magnetic phases is further moti- 
vated by applications to quantum-information processing 
and their relevance to unconventional superconductivity 
such as high-temperature and heavy-fermion supercon- 
ductivity. It has been proposed (and recently exper- 
imentally investigated [12]) that antiferromagnetic exci- 
tations may be responsible for the pairing mechanism in 
unconventional superconductors. The experimental real- 
ization of magnetic Hamiltonians with ultracold atoms 
would open up the way for direct control over inter- 
actions, geometry, and frustration. Recently an Ising 
antiferromagnetic phase has been realized with single- 
component ultracold atoms trapped in a one-dimensional 
optical lattice tilted by a magnetic field [12] (Ising density 
wave order could also be realized at lower energy scales 
in specific two-dimensional tilted geometries |14j). 

Magnetic phases with mixtures of ultracold atoms can 
be realized at low temperatures, strong interactions and 
specific filling factors. Experimentally, though, the pres- 
ence of the external confinement will challenge their re- 
alization. Previous theoretical studies, focused on Bose- 
Bose mixtures, reported accurate results for the ground- 
state phase boundaries of the Ising antiferromagnetic 
and the xy-ferromagnetic phases at half- filling factor [IS] , 
and for the critical temperature needed for their realiza- 
tion |16j . The temperature scale, determined by spin- 
exchange interactions, results challengingly low (~ few 
hundreds pK). Nonetheless, recent experimental efforts 
have been devoted towards the development of new tech- 
niques of refrigeration and thermometry [7J [5] . 

The next crucial question to theoretically address is 
studying the effect of an external trapping potential. In 
the present work we focus on the Ising antiferromag- 
netic (AF) phase realized by Bose-Bose mixtures in the 
strongly interacting regime and at half-filling factor, and 



study under which conditions it survives the presence of 
a harmonic confinement. In the local-density approxima- 
tion, the latter provides a scan over chemical potentials of 
the two components. Therefore, even if Hamiltonian pa- 
rameters are experimentally tuned to stabilize the Ising 
AF phase in the homogeneous system, the presence of 
an external confinement will challenge its realization. In 
other words, other phases such as Mott insulator (MI) + 
superfluid, double superfluid and phase separation will 
compete, moreover given the pronounced asymmetry in 
hopping parameters required for stabilization of the AF 
phase [HI EH ITS] - Hence, for a given harmonic confine- 
ment, the challenge consists in tuning the particle num- 
bers so that there exist an extended region in the trap 
with half filling for both components. 

In the following we present first precise results based on 
Quantum Monte Carlo simulations of a two-dimensional 
system with numbers of particles typical of current exper- 
iments [T7J [TS] and trapping frequencies easily realizable 
with magnetic or optical traps. We show that observa- 
tion of the AF phase in a trap is non trivial, but indeed 
possible, provided one properly tunes the total number 
of particles and the population imbalance. Hence, our 
results provide an accurate guidance towards the exper- 
imental realization of the Ising AF with mixtures. 

We consider a system of two-component bosons in a 
square lattice, with repulsive interspecies interaction and 
in the presence of an external harmonic confinement. 
This system can be realized by loading an optical lat- 
tice with two different atomic species [5], or the same 
atomic species in two different internal energy states [BJ- 
8 . If intraspecies interactions are much larger than any 
other energy scale (a limit which can be achieved by, e. 
g., using Feshbach resonance), and the temperature is 
low enough, the system is accurately described by the 
two-component hard-core Bose-Hubbard Hamiltonian: 

H = - E (*« a l a > + *«• h \ b i) - E < + u <* E 

(ij) a ;i i 

where a — a^b denotes the two species; * s intended 
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FIG. 1: (Color online) Top and middle rows: representative density maps (light and heavy component respectively) for each 
type of phase coexistence observed. Color scheme: white (black) corresponds to 1 (0) density. Bottom row: density profiles 
along a line going through the trap center (dashed red, dotted blue and dot-dashed green lines correspond to light component, 
heavy component and total density profiles, respectively). Coordinates x and y are given in units of d. Hamiltonian parameters: 
t b = 0.15t o , Uab = 5.7t a , T = 0.05t Q , w = 0.005t a . First column: 2CB core, N a = 1940, N b = 430; second column: 2CB ring, 
N a = 2155, Nb = 1030; third column: global Mott core, N a = 2130, N b = 125; fourth column: global Mott ring, iV a = 1750, 
Nb = 1750; fifth column: phase separation, N a = 425, Nb = 3185. 



on nearest-neighboring sites; aj (a]), hi (6?) are bosonic 
annihilation (creation) operators at site i, satisfying 
the hard-core constraint; t a , i& are hopping amplitudes; 
n" = a\a,i, n\ = b\bi, are number operators; U a b is the 
interspecies interaction; (if = (x a — w(ri/d) 2 , with fi a 
the chemical potential, rj the vector position of site i, d 
the lattice spacing, and w the curvature of the external 
potential. The usual intraspccics interaction term 
Hint = S Q i^a n i( n ? — 1)/^ nas been dropped because 
the U a — ^ oo limit is considered. In the following we use 
t a as the energy unit. 

In the strongly-interacting limit, Hamiltonian H can 
be mapped onto an anisotropic Heisenberg Hamiltonian 
thus realizing quantum magnetic phases [§1 HO] - In 
the following we focus on the Ising AF phase which, in 
bosonic language, is equivalent to a double checker-board 
(2CB) solid, characterized by broken Z 2 symmetry. To 
this end, we consider the hard-core boson regime set 
tb = 0.15i a , U a b = 5.7t a , so that the ground state of the 
homogenous system at half-integer filling corresponds to 
the 2CB phase [15 . We notice that the 2CB phase can 
be realized only for strong enough asymmetry between 
t a and tb, i.e. if component B is, via e.g. optical depths 
tuning, effectively much heavier than A. 

We have performed Quantum Monte Carlo simulations 



by means of the Worm Algorithm [TH] of model H for 
several curvatures of the harmonic potential w/t a = 0.1, 
0.02, 0.014, 0.007, 0.005, 0.0025, with maximum number 
of particles (corresponding to the most shallow trap) 
N = N a + N b - 7200. For each w and N the system 
size was chosen large enough so that no finite-size 
effects were present (ranging from linear size L = 60 
to L = 135, depending on the specific case), and the 
temperature was fixed to T = 0.05i a , well below the 2CB 
phase-transition temperature for the chosen Hamiltonian 
parameters [TB]. Our goal is to establish for which total 
particle number N and population imbalance (N a — Nb) 
the 2CB phase is stabilized in the trap. To this extent we 
systematically analyze and classify different typologies 
of phase coexistence. 

We have observed five distinct types of phase coexis- 
tence, as illustrated in Fig.[l] where the top (middle) row 
refers to the in trap density maps of the light (heavy) 
component, and the bottom row shows corresponding 
density profiles along a line going through the center 
of the trap. The color scheme in the top and middle 
rows is such that white (black) corresponds to density 
1 (0). In the bottom row, dashed red, dotted blue 
and dot-dashed green refer to light component, heavy 
component and total density profiles, respectively. Each 
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column is representative of an extended region in the 
plot of Fig. [2] summarizing the different types of phase 
coexistence found, as we shall discuss below. We classify 
the different typologies as follows (for each typology 
phases are listed as they appear going from the center 
to the edges of the trap): 

2CB core: 2CB, light component superfluid (Fig. [TJ 
panel al-a2). 

2CB ring: heavy component MI + light component 
superfluid, 2CB, light component superfluid (Fig. [TJ 
panel bl-b2). 

global Mott core: global MI, i.e. a MI in the total 
density (see profile of total density, green line in Fig. [TJ 
panel c3 ), light component superfluid (see Fig. [I] panel 
cl-c2). 

global Mott ring: heavy component MI + light 
component superfluid, global MI (see profile of total 
density, green line in Fig. [TJ panel d3), light component 
superfluid (see Fig. [TJ panel dl-d2). 

phase separation: heavy component MI + light 
component superfluid, heavy component MI, light 
component SF (see Fig. [TJ panel el-e2). 
For each of such sequences, subsequent phases are 
spatially separated by a thin shell of double superfluid. 
On the bottom row, density profiles make clear the 
choice of classification of typologies. In particular, in 
correspondence of the global MI phase (columns c) and 
d) ), we have observed corrugations in the density pro- 
files of the single components, which have not completely 
disappeared within simulation time. We have observed 
that the super-counter-fluid phase |10) exhibits a similar 
behavior, and suspect that the global MI phase could 
display an analogous off-diagonal-long-range-order. 

In all our simulations we have observed that the heavy 
component tends to locate itself at the center of the 
trap. This behavior can be understood as follows. In the 
strong-coupling regime considered here, derealization 
of the heavy component (at the expense of the energy 
required to climb the walls of the trap) is suppressed 
compared with the derealization of the light component, 
due to larger hopping parameter of the latter. 

Our results for the three shallower trapping potentials 
simulated (w — 0.007, 0.005, 0.0025 t a ) are summarized 
in Fig. [2j where each symbol corresponds to one of 
the typologies described above (diamond for 2CB core, 
square for 2CB ring, downward triangle for global Mott 
core, upward triangle for global Mott ring, circle for 
phase separation). We can identify (N a ,Nb) values 
(we scale the number of particles by the curvature of 
the trapping potential, see below) for which the 2CB 
phase is realized in certain regions of the trap: squares 
and diamonds, corresponding to 2CB ring and 2CB 
core, respectively. We note that, in general, the 2CB 
phase can be observed only if the light component 
is the majority one, and provided w < 0.01 t a . We 
have studied larger values of w (= 0.014, 0.02, 0.1 t a ), 
i.e. steeper traps, and have not observed CB for any 
(N a ,Nb). In other words, the trapping potential should 
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FIG. 2: (Color online) Typologies of phase coexistence (see 
text) as a function of number of particles N a , Nb scaled by 
the trap curvature for w = 0.007, 0.005, 0.0025 t a . Dia- 
monds: CB core; squares: CB ring; downward triangles : 
global Mott core; upward triangles: global Mott ring; circles: 
phase separation. Error bars are much smaller than symbol 
size (maximum percent error is 0.5%). 



be smooth enough to allow stabilization of 2CB ring or 
core regions. For w = 0.007, 0.005, 0.0025 t a , instead, 
we can summarize our results on the same plot by using 
wN a /t a and wNb/t a as axis (the smoother the trap, 
the wider 2CB core and ring regions are). Upon further 
decreasing the trap curvature w we have observed that 
2CB ring region in the wN a /t a vs wNb/t a plane tends to 
grow, although extremely slowly (e.g. lowering w by a 
factor of 2, would only convert a single upward triangle 
into a square, i. e. a global MI ring would become a CB 
ring, in the plot of Fig. [2]). In other words, for much 
shallower trap curvatures, the CB ring region would be 
only slightly larger than the one in the plot of Fig. [2] 
(i. e. the system exhibits a scalable behavior for a wide 
range of w's). Hence, Fig. 2] can be used as guidance for 
experiments since it provides N a and Nb required for the 
observation of the Ising AF for experimentally realizable 
trap curvatures. 

Although the hard-core limit corresponds to most 
stable quantum magnetic phases against hopping [9], we 
will now relax this constraint. We have performed sim- 
ulations at finite intraspecies interaction U a = Ub = U 
both in uniform and trapped system. Our goal is to 
probe how far one can soften U for given values of 
hopping parameters t a , tb, and interspecies interaction 
U a b- We have found that for tb = 0.15i a and U a b = 5.7t a 



4 



the 2CB phase in a uniform system survives down to 
U ~ 70t a , i-e. U/U a b — 12. If we further increase the 
hopping asymmetry, e.g. t b = 0.075t a , a 2CB phase still 
exists at U ~ 40i o , i.e. U/U a b — 7. We expect that 
the extension and position of 2CB regions in Fig[2] will 
slightly change upon relaxing the hard-core constraint. 
In particular, we expect the 2CB regions to shrink 
(due to a more fragile CB phase for soft-core atoms) as 
phase separation becomes more favorable. Preliminary 
numerical results agree with such an expectation but 
also suggest that the general qualitative results (i.e. 
typologies observed) valid for hard-core bosons should 
still be valid. More specifically, we have performed 
simulations in a trap with parameters w/t a = 0.005, 
t b = 0.075t Q , U ab = 5.7t a , U = 40i a , N a = 1240 and 
N b = 170 and found that the 2CB core still exist. 

A previous mean-field study carried on for soft-core 
bosonic atoms claimed realization of 2CB ring and core 
regions in a trap [5U]. The parameters used in [5D] are 
the following: w/t a — 0.064, tb = 0.l6t a , U a b = 5.7t a , 
U a = 1.9t a , Ub — 37.54- We note that component A 
is rather soft. Preliminary numerical simulations with 
the same parameters indicate that phase separation will 
occur both, with and without the harmonic trapping 
potential (in the presence of the latter complete phase 
separation, i.e. no overlap between the two components 
anywhere in the trap, seems to occur). Moreover, our 
numerical results for soft-core atoms in uniform systems, 
indicate that the 2CB phase will be destroyed upon 
softening U < 70t a with tb — 0.15t a . Hence, the results 
of [30] seem to be rather controversial. 

Finally, we would like to provide estimates for ex- 
perimental setups required for the observation of the 
Ising AF (we use a harmonic approximation around 
the minima of the optical lattice potential [3T] ). The 
hard-core limit can be achieved by realizing a a 3> a a f, 
(a = a, 6), with a Q and & a b intra- and interspecies 
scattering lengths, respectively. The latter inequalities 
can be satisfied by e.g. making use of a combination of 
dc and rf magnetic fields [22] , or dc electric fields [23] in 
order to independently tune two scattering lengths via 
Feshbach resonances. Alternatively, for what concerns 
U a b tuning, Wannier-function overlap in the presence 
of state-dependent lattices can be manipulated. Har- 
monic frequencies for the trapping potential are given 
by lj^ = 2w/(m a a 2 ). Let us consider laser beams 
with A = 1064nm and a mixture of 87 Rb atoms in 
hyperfine states 1 1 , — 1) (A component) and |2, — 2) (B 
component), for which & a b — 98.09ao (a is the Bohr 



radius), a a = 100. 4ao and a& = 98.98a . According 
to the choice in [24], with a harmonic confinement of 
uj a = u b ~ 2tt 370Hz (corresponding to w — 0.005i a ) one 
can observe the 2CB core (ring) by loading N a ~ 1700 
(2200), N b ~ 300 (800). With w a = uj b ~ 2tt 250Hz 
(corresponding to w = 0.0025i a ), instead, one can 
observe the 2CB core (ring) by loading N a ~ 3400 
(4400), Nb ~ 600 (1600). For shallower traps, e.g. 
u ~ 2tt 80Hz, one should load N a - 34000 (44000) and 
N b ~ 6000 (16000). 

In the case of 87 Rb- 41 K mixtures (we choose 87 Rb as 
the effectively-light component i. e. component A) one 
has aRbK = 163ao, aRb = 99ao and ax = 65ao. Accord- 
ing to the choice in |25j , with a harmonic confinement of 
<^Rb ~ 27r 366Hz and uk ~ 2-7T 525Hz (corresponding to 
w = 0.005i a ), the 2CB core (ring) is realized by loading 
A Rb ~ 1700 (2200), N K ~ 300 (800). With a choice of 
WRb ~ 27r 254Hz and ujk ~ 2-7T 381Hz (corresponding to 
w = 0.0025£ o ), the 2CB core (ring) is realized by loading 
A Rb - 3400 (4400), N K - 600 (1600). 

In conclusion, we have presented accurate results 
based on Quantum Monte Carlo simulations of two- 
component bosonic systems on a square lattice and 
in the presence of an external harmonic confinement. 
Starting from hopping parameters and interaction 
strengths which stabilize the Ising antifcrromagnetic in 
the homogeneous case and at half integer filling factor, 
we have studied how the presence of the harmonic 
confinement challenges the realization of such phase. We 
have considered realistic trapping frequencies and num- 
ber of particles, and established under which conditions, 
i.e. total number of particles and population imbalance, 
the antiferromagnetic phase can be observed in the trap. 
With recently developed single atom imaging techniques, 
direct observation of the Ising ferromagnetic phase is 
possible [TTl [18]. Alternatively one can use Bragg 
spectroscopy in order to probe density-density corre- 
lations [27] , Given the challengingly low-temperature 
scale associated to quantum magnetic phases, it would 
be certainly interesting to study signatures of the onset 
quantum magnetic phases on the approach of the critical 
temperature from above. 
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